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Abstract 

This paper deals with asymptotic stability of a class of dynamical systems in terms of 
smooth Lyapunov pairs. We point out that well known converse Lyapunov results for differ¬ 
ential inclusions cannot be applied to this class of dynamical systems. Following an abstract 
approach we put an assumption on the trajectories of the dynamical systems which demands 
for any trajectory the existence of a neighboring trajectory such that their difference grows 
linearly in time and distance of the starting points. Under this assumption, we prove the 
existence of a C'^-smooth Lyapunov pair. We also show that this assumption is satisfied 
by differential inclusions defined by Lipschitz continuous set-valued maps taking nonempty, 
compact and convex values. 
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1 Introduction 

Inspired by the work of Lyapunov starting in the 1950s a lot of effort has been spent on the stability 
analysis of dynamical systems in terms of Lyapunov functions. Beginning with ordinary differential 
equations defined by a continuous function, the results have been extended to differential inclusions 

x{t) G F(x(t)), x{t) G M", t > 0, x(0) = xo G M", (1) 

where F: MF R" is a set-valued map satisfying 0 G F{0). A comprehensive exploration of 
the connection between stability of differential equations and differential inclusions and Lyapunov 
functions can be found in [3] . For general references to set-valued maps and differential inclusions 
the interested reader is referred to [2] and HEU, respectively. 

A solution to ((T]) is an absolutely continuous function x : R+ —>• R” with x(0) = Xq such that H]) 
is satisfied almost everywhere. Following [71 Proposition 2.2] the equilibrium x = 0 of differential 
inclusion m is called strongly asymptotically stable if each solution can be extended to [0, oo), for 
any e > 0 there is a J > 0 such that any solution x(-) with ||x(0)|| < <5 satisfies ||x(t)|| < e for all 
t > 0, and for each individual solution x(-), one has lim 4 _,.oo x{t) = 0. 

The analysis of robust stability has been an active field in the dynamical systems literature. In 
the wake of this, the investigation of converse Lyapunov theorems and, in particular, the construc¬ 
tion of smooth Lyapunov functions is of vital interest, cf. (3(11113]. 

Clarke, Ledyaev and Stern [7] (see also [13]) have shown that, provided F{x) is nonempty, 
compact and convex for every x G R" and the set-valued map F is upper semicontinuous, i.e. for any 
X G R" and any £ > 0 there is a J > 0 such that F{y) C F{x)+eB{0, 1) for all y G x+5B{Q, 1), where 
i?(0,1) denotes the unit open ball in R", the differential inclusion ([IJ is strongly asymptotically 
stable if and only if there is a C°°-smooth and positive definite pair of functions (F, W) such that 
V is proper and 

max (VF(x),x) < —IF(x) for all x G R” \ {0}. (2) 

v&F{x) 
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A different proof of this converse Lyapunov theorem following a metric approach using weak KAM 
theory has been obtained by Siconolfi and Terrone, cf. m- Related results for retarded functional 
equations and difference inclusions can be found in [14] and [15) . respectively. 

Originating from stochastic systems, such as multiclass queueing networks and semimartingale 
reflected Brownian motions liiin], there is a class of dynamical systems that does not immedi¬ 
ately fall into framework mentioned above. More specifically, the analysis of recurrence behavior of 
the stochastic processes corresponding to multiclass queueing networks or semimartingale reflected 
Brownian motions, based on the remarkable insights of Eiiiniin], can be reduced to the stability 
analysis of a related deterministic system, called fluid network and linear Skorokhod problem, re¬ 
spectively. Both models are obtained by taking limits of scaled versions of the stochastic processes. 
In [101ITB] it is outlined that a wide class of linear Skorokhod problems and fluid networks can be 
defined by differential inclusions in a natural way. An essential part in m is to describe the linear 
Skorokhod problem in terms of a differential inclusion and construct a -Lyapunov function. The 
paper by Dupius and Williams m was published a few years before the above mentioned paper 
on smooth Lyapunov functions by Clarke, Ledyaev and Stern [7]. 

The differential inclusions considered in [7] and m are both defined by an upper semicontinuous 
set-valued map F with nonempty, compact and convex values. The techniques used to construct 
a smooth Lyapunov function have in common that the set-valued map F is embedded into a 
local Lipschitz set-valued map Fl, which keeps the property of asymptotic stability. Whereas the 
procedure in [10] uses explicitly the properties of the set-valued map describing the evolution of the 
linear Skorokhod problem, the embedding technique in [7] is applicable to any upper semicontinuous 
set-valued map taking nonempty, compact and convex values. The essential feature of the set-valued 
map Fl being local Lipschitz continuous is that it provides a Lyapunov function which is locally 
Lipschitz continuous and this property can be carried over to conclude a locally Lipschitz continuous 
Lyapunov function for the original differential inclusion. Moreover, the local Lipschitz continuity of 
the set-valued map F^ facilitates to establish that the convolution of the local Lipschitz continuous 
Lyapunov function and a C'^-smooth mollifier satisfies locally the decrease condition (|n|). The 
construction is completed by using a locally finite open covering of R” and a smooth partition of 
unity subordinate to it. 

We show that, in general the set-valued map defining the differential inclusion describing the 
evolution of a fluid network is not upper-semicontinuous. Thus, although the the zero solution 
may be strongly asymptotically stable the existence of a C°°-smooth Lyapunov pair cannot be 
concluded by the results on differential inclusions mentioned above. In this paper we follow an 
abstract point of view, starting with Zubov |26j . understanding dynamical systems as abstract 
mathematical objects with certain properties. This has been further explored by Hale, Infante, 
Slemrod and Walker, cf. [EKniEQlEl]. In the literature there are several terms used, for instance, 
generalized dynamical system, (7°-semigroup, (semi)flow, process or abstract dynamical system, cf. 
pS] and the references therein. The class of abstract dynamical systems considered in this paper is 
defined by the characteristic properties of fluid networks. The trajectories of fluid networks evolve 
in the positive orthant. In order to get a C'^-smooth Lyapunov function on the positive orthant we 
use an extension of a Lyapunov function candidate to K" by taking absolute values component-by¬ 
component. As this dehnes a continuous map the extended Lyapunov function is continuous as well. 
We note that, as the solutions to linear Skorokhod problems also stay within the positive orthant, 
Dupuis and Williams nni solved the boundary problem by shifting the orthant by some positive 
constant. Further, we note that the class of abstract dynamical systems under consideration may 
in general not be defined by a differential inclusion. As a consequence the constructions of a local 
Lipschitz continuous Lyapunov function in [adnKis], which are based on the right-hand side of 
the differential inclusion, are not applicable in the present setting. It turns out that the essential 
ingredient to obtain a local Lipschitz continuous Lyapunov function is an estimate on the evolution 
of the difference of trajectories. For this reason, we have to make an assumption on the trajectories 
of the abstract dynamical system (see assumption (A) in Theorem [T|). Considering the assumption 
from the differential inclusions perspective we show that it is automatically satisfied for every 
differential inclusion with Lipschitz continuous right-hand side. 

The paper is organized as follows. In Section 2 we state relevant notation and terminology that 
is used throughout the paper. Section 3 introduces the class of abstract dynamical systems that 
is considered and the main result is presented. In Section 4 we outline that the class of abstract 
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dynamical systems is motivated by the analysis of fluid networks. We also show that the classical 
results on smooth Lyapunov functions do not apply to the class of abstract dynamical systems 
discussed in this paper. In Section 5 we examine the relation of the assumption posed on the 
trajectories in the light of differential inclusions. Finally, Section 6 is devoted to the proof of the 
main result. 


2 Notation and terminology 

A function /: R" ^ K is called proper if the sublevel sets {x G R"|/(a::) < c} are bounded for all 
c > 0. For r > 0 and x G R" let B{x,r) := {y G R^IHx — y|| < r}. A function k G C'°°(R",R+) is 
called a mollifier if supp k = B{0, 1) and 


/ fc(x) dx = 1. 
dR" 


Furthermore, the support of a mollifier can be scaled in the following way. For r > 0 consider 
kr{x) := k{r~^x). Then, it follows kr G R+), supp kr = B{0,r), and 


/ kr{x) dx = 1. 

.n 


Moreover, to consider the convolution of a function / G and a mollifier kr, let U be 

an open subset of R” and Ur = {x £ U \d{x,dU) > r}. Then, the convolution, denoted by 
fr-Ur^ R, is defined by 


X fr{x) := f *kr (x) = / f {x - y) kr(y) dy. 

JB{0,r) 

By standard convolution results it follows fr G C°°{Ur,M.+ ), see for instance [TTl Theorem 6 
Appendix C.4]. Furthermore, if / is continuous in U, it holds fr ^ f uniformly on compact 
subsets (u.o.c.) of 17 as r —>• 0. The Dini subderivative of a function / :17—^-RatxGl/in the 
direction v G R" is defined by 

C/(x;»):= liminf "/M . 

5 —£ 


Let T(x,R") denote the contingent cone to R" at x defined by 


r(x,R!;:) = <{ n G R” I liminf + ^ g 

^ I £->-0 £ 


with d{x, K) = inf{||x — y\\ \y G K}. 


3 Statement of the main result 

We start this section by recalling an abstract definition of a dynamical system from [^. A 
dynamical system defined on a metric space X is a continuous mapping u: R+ x A —> A such that 
u(0, x) = X and 

u{t, u{s, x)) = u{t + s, x) for all t,s € R+, x G A. 

Recall that x* G A is an equilibrium if u{t,Xr) = x* for alH > 0 and x* is said to be stable if for 
every e > 0 there is a i5 > 0 such that d{x, x*) < 6 implies that d {u{t, x), u{t, x*)) < e for all t > 0. 
If in addition, there is a M > 0 so that d(x,x*) < M implies that hmt_^oo d{u{t,x),u{t,x^)) = 0, 
then is called asymptotically stable. 

Here we consider the metric space V C C(R+,R") defined by the following properties: 

(a) There is a L > 0 such that 

||:/5(t) — V?(s)|| < L\t — s\ for all ip € V, t, s £ R+. 
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(6) Scaling invariance: ^ (p{r ■) gV for all (f GV,r > 0. 
(c) Shift invariance: ip{- +t) gV for a\\ ip GV,t > 0. 


(d) If a sequence (v5n)neN in 'P converges to uniformly on compact sets, then ip^ G V. 


(e) 


Concatenation property: 
Ot* ip 2 where 


For all ipi,ip 2 G P with ipi{t*) = 952 ( 0 ) 


for some t* > 0 it holds 


ipi Of ip2{t) := 


Vi{t) 

Mt-t*) 


t < t*, 
t > t*. 


(/) There is a T > 0 such that the set-valued map P : M" V defined by 

P{x) = {^ : [0, T] ^ R” I y: e V, <^(0) = x} 


is lower semicontinuous, i.e. for any ip G P{x) and for any sequence of elements (x„)„gN G R+ 
converging to x, there exists a sequence (99n)rtgN with ipn G P(xn) converging to ip uniformly 
on compact sets. 


By condition (a) the functions ip G V are Lipschitz continuous with respect to a global Lip- 
schitz constant. Condition (c) is in one-to-one correspondence to time-invariance of differential 
equations/inclusions. Condition (d) expresses that the set V is closed in the topology of uniform 
convergence on compact sets. The concatenation property (e) and condition (f) imply that for 
X G R" the set of functions in V starting in x depend lower semicontinuously on x, cf. m- 

In order to define the class of dynamical systems which will be considered in this paper, we 
equip the set P with the metric 


d(iPi,iP 2 ) ■= max 2 
NeN 


-N 


Il9?i - 9^21^ 

1 -I- II 951 - (/ 92 IU’ 


where ||99||jv := supjgjQ^^j ||v5(t)|| so that convergence of functions is equivalent to uniform conver¬ 
gence of the corresponding restrictions on each compact subset of R+, cf. m- Furthermore, we 
define the shift operator 


S{t): C(R+,R!^) ^C(R+,R”), S{t)ip{-) = +1). 

In this paper, we consider the dynamical system on the metric space P C C'(R+,R") defined by 
the mapping 


u: R+ xP ^ P, u{t,ip) = S{t)ip{-) = ip{- -I- 1). 

Throughout the paper we call a function ip G P a trajectory of P. The zero trajectory 95 * = 0 
is the unique fixed point of the shift operator S{t) defined on P and thus, 99 * = 0 is the only 
equilibrium of the dynamical system. The scope of this paper is to characterize asymptotic stability 
of the dynamical system u defined on P in terms of the existence of a smooth Lyapunov function. 

A pair {V, W) of positive definite functions on R” is called a Lyapunov pair for the dynamical 
system u defined oaP liV ■. R” —>■ R+ is proper and for any ip gP, 

V{ip{t)) — V{ip{s)) < ~ / ff^(95(?')) dr for all 0 < s < t G R+. (3) 

J S 

The pair [V, W) is called a C°°-smooth Lyapunov pair if the functions V and W are C^-smooth. 
To formulate the decrease condition in a differential form let 

+ (4) 

The main result of the paper is the following. 

Theorem 1. Suppose the dynamical system u defined on P satisfies: 
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(A) For any ip € P, e > 0, and T > 0 there is a continuous function c: [0,r] —>■ ]R_|_ such that 
limt^o exists and is positive and for any y £ K." with y)(0) — y £ B{p{0), e) flR" there is 
a trajectory if £V with ip{0) = p{0) — y satisfying 

II Fit) -y- tPit) II < ||y|| c{t) for all t £ [0, T], 

Then, the dynamical system u defined on V is asymptotically stable if and only if there is a C°°- 
smooth Lyapunov pair (V, W) such that for every p £V it holds 

^(Fit)) < ~Wipit)) for all t>0. (5) 

We note that estimates similar to (A) for trajectories of differential inclusions with state con¬ 
straints were derived by Bressan, Facchi, Bettiol, Vinter and Rampazzo, cf. [5] and the references 
therein. 


4 Motivation and application of the main result 


For the purpose of this paper we give a very brief introduction to fluid networks. For a com¬ 
prehensive description of multiclass queueing networks and fluid networks we refer to H 0 ig. 
A fluid networks consists of J € N stations serving n € N different types of fluids with J < n 
and each fluid type is served exclusively at a predefined station. This assignment defines the 
constituency matrix C £ with Cjk '■= 1 if fluid type k £ n} is served at station 

j £ {1,..., J} and Cjk = 0 otherwise. The exogenous inflow rate of fluid type k is denoted by au and 
a = (ai,..., Q!„)^ £ R" := {x £ : Xi > 0, i = 1,..., n} is called the exogenous inflow rate. Like¬ 

wise, pk G R+ denotes the potential outflow rate of type k fluids and p = {pi, ...,pnj^ £ R+. The 
substochastic matrix P £ [0, describes the transitions in the network, where it is assumed that 
the spectral radius of the matrix P is strictly less than one, i.e. 1 > max{|A| | 3 a: ^ 0 : Px = Aa;}. 
The initial fluid level and the fluid level at time t of the network are denoted by a;(0) and x{t), 
respectively. We note that, as x{t) described the deterministic analog of the queue length of the 
multiclass queueing network, the fluid level process a;(-) evolves only in the positive orthant R”. 

The discipline, for instance First-In-First-Out (FIFO), determines the rule under which the 
individual stations of the fluid network are serving the different fluid types. More formally, the 
discipline specifies the allocation rate, denoted by u = (ui,...,u„)^ £ R", where Uk denotes the 
current amount of time that a predefined station allocates to serve fluids of type k. For the class 
of general work-conserving fluid networks, given x £ R" , the set of admissible allocation rates is 

C/(a;) = {it G R" I M > 0, e — Cu > 0, (Cx)~^ ■ (e — Cu) = O}, (6) 


where e = (I,..., 1)^ £ R*^ and the inequalities have to be understood component-by-component. 
The evolution of the fluid network can then be described by the following differential inclusion 

X £ G(x) := |a — (I — P^)Mu | u £ fl T(a;,R"), (7) 

where M = diag(/i) and r(a:,R") denotes the contingent cone to R" at x. The intersection with 
the contingent cone to the positive orthant is to ensure the nonnegativity of the solutions. Tackling 
a simple example we show 

Lemma 1. The set-valued map G defined in © is not upper semicontinuous in general. 

Proof. To show the claim we consider a single station fluid network serving one type of fluid. That 
is, for a > 0, ^ = a 3- 1, and P — 0 the differential inclusion © is defined by set-valued map 

G(a:) = {a — (a-I-I) 11 I 0 < it < 1, a:(I — it) = 0} fl r(a:, R+) 

f I I 

= {q: — (a-l-l)it I 0<u<l,it=lifa;>0 and it = else > . 

To conclude that G is not upper semicontinuous let a; = 0 and consider the sequence ixk)ken with 
Xk = ]:■ Then, for each k £N we have Gixu) = — I. Let graph(G) := { {x, y) G R” xR” | y G G(a:) } 
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denote the graph of G and consider the sequence (xfc, — on the graph of G which converges 

to (O, — l) ^ graph(G'). Hence, the graph is not closed and by Proposition 2 in [J Section 1.1] the 
set-valued map G is not upper semicontinuous. □ 

We note that the differential inclusion defined by ([8|) is strongly asymptotically stable. However, 
as a consequence of Lemma [H the existence of a C°°-smooth Lyapunov pair cannot be concluded 
by the results on differential inclusions mentioned above. 

We use the main result Theorem [T] to conclude that the differential inclusion (|8]) admits a G°°- 
smooth Lyapunov pair. It is well known that a single station general work-conserving fluid network 
serving only one fluid type satisfies the properties (a)-(f); cf. [H1[TH1H5]. To conclude the existence 
of a C'^-smooth Lyapunov function we will apply the main result of the paper by interpreting the 
zero solution </?* = 0 of the differential inclusion © as the fixed point of the shift operator defined 
on the set of solutions to 0. 

Theorem 2. The differential inclusion (jS)) admits a C°°-smooth Lyapunov pair, i.e. there is a 
C°°-smooth and positive definite pair (V, W) such that V is proper and 

max {W(x),v) < —W(x) for all x G (0,oo). (9) 

«eG(x) 


Proof. To show the existence of a (7°“-smooth Lyapunov pair we verify that the set of solutions 
Sg to the differential inclusion ([HI satisfies the assumption (A). 

Let ip G Sg, e > 0 and T > 0 be fixed. In a first step, we treat the case that (/3(0) > 0. Then, 
we have 


p{t) 


(p(0) — t if t < (p(0) 
0 else. 


In the case y = (/3(0) the only solution if to the differential inclusion starting in ■0(0) = (/3(0) —y = 0 
is the zero solution and we obtain \p{t) —y — tp{t)\ = jtj for all t < (^(0) and \p{t) —y — tp{t)\ = \y\ 
otherwise. Hence, one has 

\p{t)-y-fi{t)\<^^\y\t foralH>0. 

If ?/ ‘/^(O) ’"^6 consider the solution ip to the differential inclusion starting in ^’(0) = </5(0) — y given 

by 




‘^(0) -y-t if t < (^(0) - y 
0 else. 


On one hand, if 0 < <^(0) — y < (/?(0) it follows \p{t) — y — ip{t)\ = 0 for all t G [0, (/?(0) — y\. Also, 
for all t G [‘/?(0) — y, ‘/^(O)] we have \p{f) — y — ip{t)\ = |‘/?(0) — y — t\ and for all t > (/?(0) one has 
~ y ~ '0(^)1 = |j/l- On other hand, if 0 < (/^(O) < v?(0) — y we have \p{t) —y — 'f{t)\ = 0 for all 
t G [0,1^(0)]. Further, for all t G [v5(0), i^(0) — y] one has \p{t) — y — '0(t)| = ||^(0) — y — t\ and for 
all t > 1^(0) — 2 / it follows \p{t) — y — fj{t)\ = \y\. Consequently, in either case one obtains 

\pit)-y-f’{t)\<^^\y\t foralH>0. 

Finally, we consider the case (/^(O) = 0. Then, we have = 0. For y G (—e, 0] a solution f} of the 
differential inclusion with tffiS) = —y is 




-y-t \it<-y 
0 else. 


Therefore, \p{t) —y — V'(0l = ^ ^ ^ Pj “2/1 \p{t) —y — tp{t)\ = \y\ for all t > —y and one 

has 

—y — fi’it)\ < log(t -I- e) \y\ for all t > 0. 
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Thus, there is a C'°°-smooth Lyapunov pair {V.W) such that for every solution ip € Sc one has 


Vipit)) < —W{p{t)) for all t>0. 
Consequently, the pair {V,W) satisfies 


max {V{x),v) < —W{x) for all x G (0, oo). 

vGG(x) 


This shows the assertion. 


□ 


5 Relation to differential inclusions 

Due to the fact that Clarke, Ledyaev and Stern [7] as well as Dupius and Williams [in] embed the 
set-valued map defining the differential inclusion into a Lipschitz one, we consider the differential 
inclusion 


x{t) e F{x{t)), (10) 

where F : R” R” is Lipschitz continuous, i.e. there is a constant L > 0 such that 

F{x) C F{y) -I- L ||x — y\\ i?(0,1) for all x,y G R", 

and F{x) is nonempty, compact and convex for every x G R". Let Sf{x) denote the set of solutions 

p: R+ —)• R” to (IT^ with (/3(0) = x. Let Sp ■= | x € R”}. Obviously, the mapping 

u : R+ X Sf —j' Sf, u{t, p) := p{t + •) (11) 

defines a dynamical system on 5^- 

Next we show that condition (A) is a natural assumption which is satisfied by a wide class of 
dynamical systems. To this end, we investigate assumption (A) from differential inclusions per¬ 
spective and show that it is automatically fulfilled for differential inclusions defined by a Lipschitz 
continuous set-valued map. 

Theorem 3. Let F be a Lipschitz continuous set-valued map taking nonempty, compact and convex 
values withO G F(0}. Then, the dynamical system u defined onSp by (fTTl) satisfies assumption (A). 

Proof. Let p G Sf, £ > 0 and T > 0. We define c{t) := e^* —1. Then, for y e R" the function 
Py{-) '■= p{-) — y is absolutely continuous with Py{0) = 93 ( 0 ) —y. Further, as F is Lipschitz it holds 
F{py(t)) C F{p{t)) L ||y|| B{0, 1) and we have 

d{py{t), F{py{t)) =d{p{t),F{py{t)) <L||y||. 

Thus, by Filippov’s Theorem [TJ Theorem 1 in Chapter 2, Section 4] there is a solution to (1717)1 
defined on the interval [0,r] with tj}{0) = Py{0) = (p(0) — y satisfying 

\\py{t) - i;{t)\\ = \\p{t) - y - 'fit)]] < \\y\\ (e^‘-l) for all t G [0,r]. 

This shows the assertion. □ 


6 Proof of the main result 

Before proving Theorem [T] we present a useful characterization of asymptotic stability for the class 
of dynamical systems under consideration which is based on the scaling property (b) and the shift 
property (c). 

Lemma 2. The dynamical system u defined on V is asymptotically stable if and only if there is a 
T > 0 such that u(||i^(0)||t, 1 ^) = 0 for all p GV. 
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Proof. Suppose there is a r < oo such that M(||(/3(0)||r, (/?) = 0 for all (f gV. To conclude stability 
let e > 0 and 6 := where for cc € K let [cc] := min{fc G Z | fc > a;}. Let ip gV with d{0, p) < 5. 
By assumption, using the scaling and shift property it holds 

p{s + t) = 0 for alH > r ||(/3(s)||. (12) 

Together with the Lipschitz continuity of p for every t G [0, t||(^(s)||] we have 

+ t)\\ = II p{s + t)- p{s + t||v3(s)||) II < L |t - t||(^(s)|| | < Lt||(^(s)||. (13) 

By (HH) we conclude that m holds for all f > 0. Therefore, 

ll‘^(• + 0IU= sup ||v5(s + t)|| < rTT]||(^||Ar = |||■Lr]v5||JV. 
se[o,Af] 

In turn, by the triangular inequality, we have 

d{0Mt,<p)) < = d(0, \Lr]p) < \Lr]d{0,p) < e. 

By assumption, it holds v5(||y)(0)||T+t) =0forallf>0. This in turn implies lim(_>oo d(0, u(t, (^)) = 
0 and we have attractivity. 

Conversely, let = 0 be asymptotically stable. Due to the scaling property it suffices to 
consider trajectories p with ||</?(0)|| = 1. Then, as 

lim d(0, u{t, p)) = lim d(0, p{t + •)) = 0 

t—^oo t^OO 


we have 

lim ||</?(f)|| = 0 for all p gV. 

t—^OO 

Hence, inf{||(^(f)|| 11 > 0} = 0 for any p G V with ||</5(0)|| = 1. The assertion then follows from 
Theorem 6.4 in [^. □ 

Proof of Theorem [T] In [181 Theoerm 2] it is shown that for the dynamical system under con¬ 
sideration there is a continuous Lyapunov pair if and only if the dynamical system is asymptotically 
stable. Thus, the non-converse implication is already shown. 

Conversely, let the dynamical system u defined on V be asymptotically stable. Then, by 
Theorem 2 in [H] there is a continuous Lyapunov pair (V, W) such that 

V{p{t)) — H((^(s)) < — / W{p{r)) dr for d\\ p G V, 0 < s < t. (14) 

J S 

Thus, the construction of a C°“-smooth Lyapunov-pair remains. 

To get differentiability on the boundary of the positive orthant, we first extend the pair {V, W) 
to K”. To this end, let | • |vec denote the map that takes componentwise absolute values defined by 
|a;|vec := (kiL kn|)^ G R+- The extention of the pair (H, W) to R" is dehned by 

V%x) := H(|cr|vec), W<^{x) := W{\xU)- 

Note that, as a composition of continuous functions, the pair (I4®,IT®) is also continuous. As a 
Hrst consequence of assumption (A) we conclude that is locally Lipschitz. 

Lemma 3. Suppose the abstract dynamical system u defined on V satisfies (A) and is asymptoti¬ 
cally stable. Then, is locally Lipschitz on R". 

Proof. Let [/ C R" be open, convex, and bounded and let x G U. Following Corollary 3.7 in [S], 
since is lower semicontinuous, it suffices to show that there is a M >0 such that for any 

u G R" it holds 


D{-V^){x]v) < M||u||. 



Let v' G M" and ^ > 0. Let (f gV he a, trajectory of the dynamical system satisfying <p(0) = |a;|vec 
and 

r°o rW^W-r 

V‘^{x)= \\ip{s)\\ds= ||</5(s)||ds, 

Jo Jo 

where in the last equality used the stability of V and Lemma[^ By the continuity of | • |vec we have 

lim |vec — 1^1 vec- 

C ->-0 

So, for every e > 0 and ^ sufficiently small we have |x + ^r)'|vec G i3(|x|vec,£) HR” . Moreover, there 
is a continuous mapping g : K" —>• K” satisfying ||g(r^')ll = Ik^ll 

\x |vec — l^lvec “t“ )■ 

For T := max{||a;|| r, ||a: + r} < T(||a;|| + e), by assumption (A) and the triangular inequality, 
there are c > 0 and ^/> G 7^ with ?/;(0) = |a;|vec + ^9{v') such that 

II ^(t) II - II m II < II m - m II < ciiffK)ii (1+ c{t)) = i iiu'ii (1+ c[t)) (15) 

for all t G [0,T]. The definition of the stability of V together with Lemma [U and || |a:|vec + 

^ffK)ll = l|a; + ^^^'ll yield 

/■oo /■lk+{«'lh 

V^'ix + ^v')^ ||^(s)||ds=/ ||V’(s)||ds. 

Jo Jo 

On the one hand, if ||x|| < ||a; + ^u'|| by using (ITBI) it follows 


V^ix) - V^ix + ^v') < f II ip{s) II ds - f 
Jo Jo 


Ik+C-u'lh 


II V'(s) II ds 

< / Il75(s)|| - IIV’(s) II ds 

Jo 

< [ ^ Ik'll (1 + c(s))ds < elk'll • Ikll rC, 

Jo 

where C := max (1 + c(s)). On the other hand, to consider the case ||a;|| > ||x + we note 

0<s<||x||t 

that the triangle inequality together with the Lipschitz condition imply 


||V5(i)|| < 117^(0)11 +Lt< ||x|| (1 + Lt) for all t G [0, ||x||t]. 
Using (fTSl) . (ITCl) . and 0 < ||x|| — ||x + < ilk'll we obtain 


(16) 


V^{x) - V^{x + iv') < f II (^(s) II ds - f 
Jo Jo 




llk(s) II ds 


< 


IU+5w'||r 


II 75 ( 5)11 - llk(s) I|ds + 


' IU+5«'|h 


II 75 ( 5 ) II ds 


q|x+{t;'||T 


< 


e Ik'll (1 + c(s)) ds + r( ||x|| - ||x + e^^'ll ) • sup || (p{s) || 

se[||a;+{t;'|U.IUI|T-] 

< elk'll \\x + iv'\\TC + re Ik'll Ikll (1 + Lt). 


Consequently, taking limits and using that U is bounded there is a M > 0 such that 

Di-V‘^)(x;v) = liminf M ^ ^ < r (C + 1 + Lr) ||x|| ■ ||u|| < M ||r>||. 

>11 e 

The shows the assertion. □ 
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Proceeding with the construction of a smooth Lyapunov pair, let U be an open subset of M" 
and consider the convolution of 1^® and kr defined by 

V^{x) :=V^ *kr{x) = [ V^{x-y)kr{y) dy = [ V{\x - y\^ec) kr{y) dy. 

Jr" Jr" 

Also, we consider the convolution of bP® and kr given by 

W^{x) := bP® * kr (cc) = [ W%x-y) kriy) dy. 

JR" 

By standard convolution results it follows IP® G C°“(C/, R+) and bP® G C°°(17, M+). Furthermore, 
since P® is continuous on U it holds 1/^® —)• P® uniform on compact subsets of C/ as r —>■ 0. 
Consequently, for every e > 0 there is an tq such that for all r G (0,ro) we have P/ and bP® are 
smooth on U and 

|P/(a;)-P®(a;)| <e, |bP;(a;) - bP®(a;)| < § for all a: G t/. (17) 

The subsequent statement addresses the decrease condition of the convolution along trajectories 
ip GV. 

Lemma 4. Let U C be compact such that U Cl M" ^ 0 and suppose (P, bP) satisfy (fTH) and 
assumption (A) is satisfied. Then, for every e > 0 there exists a ro > 0 such that for all r G (0,ro) 
we have 


P;(^(t))<-bP®(<p(t))+e (18) 

for all ip and t G [0,r] with ‘/3(-)l[o,T] C U flK". 

Proof. Let tp G "P be a trajectory satisfying (/5(0) = a; G 17 fl M". Then, for /i > 0 we have 

Vf{ip{t + h)) - Vf{ip{t)) = f ( V^{ip{t + h)-y)- V^{ip{f) - y) \ kr{y) dy. 

There is a continuous mapping ^ —)■ satisfying || 5 (?/)|| = \\y\\ and 

\T{t) - ylvec = ip{t)-g{y). 

Further, by assumption (A) and for h sufficiently small there is a continuous function c: K+ —^ K+ 
with lim(_>.o ^ =; cq > 0 and a trajectory '!/;(t + •) G P with 'if^t) = ip{t) — g{y) such that 

\\ip{t + h) - y - tlj{t +h)\\ < \\y\\c{h). (19) 

Using this, as V^{tp{t + h)) = V{if{t + h)) we obtain 


Vf{ip{t + h)- Vf{ip{t) < [ P®(<p(t + h)-y) - P®(V'(t + h) kriy) dy 

+ [ (viipit + h)-Vi(pit)-g{y))'\kriy) dy. (20) 

By the local Lipschitz continuity of V with constant L and (113, the first term on the right hand 
side in the above inequality can be estimated as follows 


P®((/?(t+ /i)-y) - P®('i/’(t + /i)) kriy) dy 

= f ViTit + h)-giy)) - Viipit + h)) kriy) dy 
JR^ 

< / L\\ipit + h) - giy) - ifit + h)\\kriy) dy 
JR-^ 

< cih) L f \\giy)\\ ■ kriy) dy = c(/i) L f 
JR-^ JR^ 


kriy) dy- 
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Furthermore, it holds 
follows 


'B(O.r) 


\\y\\kr{y)dy < 


'B(O.r) 


rkriy) Ay 


r and choosing rp 


£ 

2 cq L 


it 


f + h) -y) - V{i}{t + h))\ -kriy) d 2 /< 

JR^ 


Asymptotic stability of the dynamical system implies that the very last term in (uni) can be esti¬ 
mated by means of the function W and its mollification, 




< 




( V(ip(t + h)) -V((p(t) - g{y)) ) kr{y) dy 

\ 

W{iIj{s)) ds ) • kr{y) dy 

W{jp{t +s)) ■ kr{y) dy'] ds, 

/ 

where the last identity is obtained by integration by substitution. Next, we show that the function 

[ W{ip{t +s)) kr-{y) dy 


is continuous in [0, h]. To see this, consider the modulus of continuity of the function 

s i-5> W{if{t + s)), 


defined for 5 € [0, h] by 


m((5, VF(V’(t + ■))) 


sup 

|s-s'|<5 


W {f}{t -b s)) - IF [if^t + s ')) 


Then, for s, s' G [0, h] it holds 

W(ip{t + s))) -W[f}{t + s')) < m^h, W[f}(t + •)))• 

By asymptotic stability of the dynamical system Wipit + s)|| is bounded and, hence, W{'tp{t + •)) is 
uniformly continuous. Thus, we have 

lim m w((ti){t -b = 0. 

That is, for every e' > 0 there is a 5' > 0 such that m^/i, W ( 'ip{t + •))) ^ ^ h < 6'. For 

e' > 0 choose <5 > 0 such that |s — s'| < ^ < 6'. Then, 


-bs)) -W{'if{t +s'))] kr{y) dy 


< [ m(6,W{'ip{t +■))] kriy) dy < [ s'kr{y) dy = e'. 
Moreover, by conditions (El) we have —Wf{x) + § < — VF®(x) -b e. Finally, the collection of the 
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above relations yields 


h—^0 h 

-1 “ i L [l 

- ~ 1-0 1L (L ”” J 

= -[ W{'ilj{t))kr{y) dy+^ = - [ W{(p{t)-g{y))kriy) dy+^ 

= -[ W^(p{t) - y) kr{y) dy + ^ 

JR'^ ^ 

= -W,-{^{t)) + ^<-W^^{t))+e. 


This shows Lemma U] □ 

Now, let lA = {Ui}^i be a locally finite open cover of R” such that for every i the closure Ui 
is compact. Further, let {ipi}'^! be a smooth partition of unity that is subordinate to U. Define 

£i = 7 min{min y®(a:), min r(;®(a:)} and = max || V'(/'i(x)||. (2i) 

x^Ui x^Ui x^Ui 

Then, by Lemma|l]for every i there is a (7°“-pair {Vf, Wf) such that for every x £ Ui, 

\V‘^{x)-V:^{x)\ < 2.+i(iVg,) |fT^(cr)-lF,^(x)| <£,. (22) 


Moreover, by the conditions (IT^ and m we have that 

(23) 


Next, we define 


Vsix) ■.= ^il}^{x)V,^{x). 
2=1 


The following estimate holds true 


\y!{x) - V‘^{x)\ < ^ V’i(a^) \Vi{x) 

i=l 


r"(x)| < 


V<^{x) 


E 

2 = 1 




2*+i(l + g,) 




Using the triangular inequality, the latter estimate shows that U/ is proper and positive definite. 
The next step is to derive that U/ is decaying along trajectories of V. To this end, we consider 


^ [v^y^it))+ u;(v.(t)) - 


= air(^W)] + al 


oo 

.i^l 


CO CO 


2=1 


2=1 

CO 


i=i 
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Using the conditions and (HI) we get the following estimate 




oo 



Defining Si := max{ej : x € Ui (lUj ^ we have that 



OO 


= +ei). 


Using m and the triangular inequality applied to the second inequality in (l22)l . it holds that 




Finally, we have that 


OO 




Consequently, the pair {Vg,Wg) defines a C°°-smooth Lyapunov pair, which shows the asser¬ 


□ 


tion. 
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